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Abstract 

In this paper, we study the growth index of matter density perturbations for the power law 
model in f{T) gravity. Using the parametrization 7(2;) = 7o + 7i for the growth index, which 
approximates the real evolution of j{z) very well, and the observational data of the growth factor, 
we find that, at la confidence level, the power law model in f{T) gravity is consistent with the 
observations, since the obtained theoretical values of 70 and 71 are in the allowed region. 
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I. INTRODUCTION 



Cosmological data from a wide range of sources have indicated that our Universe is 
undergoing an accelerating expansion l|-|3|]. Basically, two kinds of alternative explana- 
tions have been proposed for this unexpected observational phenomenon. One is the dark 
energy with a sufficient negative pressure, which induces a late-time accelerating cosmic 
expansion. Currently, there are many candidates of dark energy, such as the cosmological 
constant, quintessence, phantom, quintom, and so on. The other is the modified grav- 
ity, which originates from the idea that the general relativity is incorrect in the cosmic 
scale and therefore needs to be modified. Examples of such theories are the scalar-tensor 
theory, the f{R) theory and the Dvali-Gabadadze-Porrati (DGP) braneworld scenario, et 
al.. 

Recently, a new interesting modified gravity by extending the teleparallel theory 4|, 
called /(T) gravity, is proposed to explain the present accelerating cosmic expansion 5|. 
Since the teleparallel theory is based upon the Weitzenbock connection rather than the 
Levi-Civita one, it has no curvature but only torsion. In analogy to the well-known f{R) 
gravity obtained from extending the Einstein-Hilbert action to be an arbitrary function 
of R, the fjT) theory is built by generalizing the action of the teleparallel gravity to be 
T + /(T) 5|, l6|. An important advantage of the /(T) gravity is that its field equations 
are second order as opposed to the fourth order equations of the f{R) gravity. So, it has 
spurred an increasing deal of interest in the literatures [7-151 . For example, some concrete 

n n n n 

models are built, in [5|, l6|, llO|, to account for the present cosmic expansion and the 
models with the phantom divide line crossing are proposed in 10, lUl]- The observational 
constraints on model parameters are discussed in 12] and the dynamical analysis for a 
general f{T) theory is performed in jisl . 

It has been pointed out [g], [i2| that the /(T) theory can give the same background 
evolution as other models, such as the ACDM and DGP. So, in order to discriminate the 
/(T) gravity from other models, one needs to break the degeneracy of the background 
expansion history. An interesting approach to differentiating the modified gravity and 
dark energy is to use the growth function 6{z) = Spm/pm of the linear matter density 
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contrast 
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361]. While different models give the same late time expansion, they may 



produce different growth of matter perturbations 23j. To discriminate different models 
with the matter perturbation, usually, the growth factor g = is used and it can be 
parameterized as 0, 



25| 



9 = 



dliaS 
din a 



(1) 



where 7 is the growth index and VLm is the fractional energy density of matter. This 



approach has been explored in many works 
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-361, and it is found that different models 



may give different values of the growth index if a constant 7 is considered. For example 



7. 



0.5454 
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for the ACDM model and 700 ~ 0.6875 [26|, \2J\ for the flat DGP 
model. Therefore, in principle, one can distinguish the modify gravity from dark energy 
with the observational data on the growth factor. 



Recently, the authors in Refs. [l4 



15| have discussed the matter perturbations in 



the /(T) gravity and found that, despite the completely indistinguishable background 
behavior, the growth of matter density perturbation of different models can be different. 
In this paper, we plan to investigate the growth index of matter density perturbations 
in the f{T) gravity, and then by using a feasible parametrization of the growth index to 
compare the growth function with the observational data. Our results show that the /(T) 
gravity is consistent with the observations. 



II. THE /(T) THEORY 

In this section, we give a brief review of the f{T) gravity. The torsion scalar T in the 
action of the teleparallel gravity is defined as 

T ^ S^T^^, , (2) 

where 

Sr = \ {K^\ + ^^T"'; - 5:T'^^^) , (3) 



and T"^^^ is the torsion tensor 



= ^A{d,et - d^et) ■ (4) 
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Here is the orthonormal tetrad component, where A is an index running over 0, 1, 2, 3 
for the tangent space of the manifold, while n, also running over 0, 1, 2, 3, is the coordinate 
index on the manifold. K'^'^^ is the contorsion tensor given by 

(5) 

Similar to f{R) gravity, the /(T) theory is obtained by extending the action of teleparallel 
gravity to be T + /(T) to explain the late time accelerating cosmic expansion with no 
need of an exotic dark energy. 

Assuming a flat homogeneous and isotropic Friedmann- Robertson- Walker universe de- 
scribed by the metric 

ds'^ = di^ — a^{t)5ijdx^dx^ , (6) 

where a is the scale factor, one has, from Eq. ([2]), 

T = -6H^ , (7) 

with H = a/a being the Hubble parameter. Thus, in f{T) gravity, the cosmic background 
equation can be written as 

H'~P-{- ^H^fT , (8) 
._ 1 6H' + f + l2H'fT 

4 1 + fT- fTT ' ^ ' 

with J't = df / dT. Here we assume that the energy component in our universe is only the 
matter with radiation neglected. Apparently, the last two terms in the right hand side 
of Eq. (IH]) can be regarded as an effective dark energy. Then, its effective energy density 
and equation of state can be expressed, respectively, as 

Pe,,-^(-/ + 2T/.) (10) 

f/T-fT + 2TfTT . . 

"^'^ (l + /T + 2T/rT)(//T-2/T) ■ ^ ' 
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In order to exp 
proposed in Refs 



ain the present accelerating cosmic expansion, some f{T) mode 
11| . In this paper, we only consider a power law model 



3, y, Q, 



/(T) = a(-T)" = a{6H 



2\n 



s are 



,l6| 



(12) 



with 



a = ml?-"' 



1 — f2mo 
2n-l ' 



(13) 



where VL^ 



mO 



87rGp(0) 



is the dimensionless matter density parameter today. Substituting 



above two expressions into Eq. (8) and defining E'^ = H^/Hq, one has 



= Qrnoil + Zf + (1 - nmo)E 



2n 



(14) 



The reason to consider the power model in our paper is that it has the same background 
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53l ] and it reduces to the ACDM 
1/2. In addition, it has a 



evolution equation as some phenomenological models 
model when n = 0, and to the DGP model 5J] when n 
smaller xitin value than the ACDM when fitting the recent observations such as the Type 
la Supernova (Sne la), the baryonic acoustic oscillation (BAO) and the Cosmic Microwave 
Background (CMB) radiation [l^. Let us note that, in order to be consistent with the 
present observational results, it is required that |n| ^ 1 



III. GROWTH INDEX OF /(T) MODEL 

To the linear order of matter density perturbations, the growth function 6{z) at scale 



much smaller than the Hubble radius satisfies the following equation 



24| 



S + 2H6-4:7rGeffpm5 = 0, (15) 

where Gefj is the effective Newton's constant and the dot denotes the derivative with 
respect to time t. In the general relativity, Gg// = Gn, where Gn is the Newton's 
constant. Defining the growth factor g = dln6/dlna, one can obtain 
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In general, the analytical solution of above equation is very difficult to obtain, and thus, 
we need to resort to the numerical methods. 



In the f(T) gravity, Gg// can be expressed as [15 1 

So, the growth factor satisfies the following equation: 

H .\ 3 



(17) 



1 



Using 

and Eq. f lT5]) . we get 
dg 



din a 



+ 9^ + 9 



H 

3 1-E 



3 1 + f/6H^ + 2fT 

'21 + fT-l2H^fTT' 



2ri-2/ 



1-Q 



mO) 



where Qr 



21-nE^^-^{l-n^o) 

z)^^ and E satisfies 
dE^ -3^2 + 3E2'^(1 - n^o) 



2l 



n(l-Q„o)g2"-2 



(19) 



(20) 



2n-l 



(21) 



din a 1 -nE2"-2(l -r]„o) 

Thus, given the values of fimo and ra, the value of g can be obtained by solving Eqs. fl20|) 
and (!2T]) numerically with the initial condition g = 1 aX z oc. Then, using Eq. ([1]), 
we can get the value of growth index. The results are shown in Fig. ([T]). From the 
left panel of this figure, we can see that, when n = 0, our result reduces to that in the 
ACDM. As expected, the growth factor g in f{T) gravity {n ^ 0) grows slower than 
that in general relativity, which is the same as that obtained in Refs. jl^ \\\ . because a 
weak effective Newton gravity (see Eq. (fT7|) ) is obtained. The right panel shows that the 
growth index 7 in /(T) gravity is different from that in the ACDM at both low redshifts 
and high redshifts. This is caused by the fact that the /(T) gravity and ACDM give 
different background evolutions at low redshifts. However, 7 is mainly determined by the 
low redshift evolution, since, at high redshifts fi^ ~ 1? different values of 7 may give the 
same result for g. When fi^o = 0.272, the values of growth index 7 approach to 0.58 and 
0.71 with z ^ 00 for n = 0.1 and 0.25, respectively. These results are different from that 
obtained from the ACDM and DGP models where the values are 700 ~ 0.5454 26|, |28| and 
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7oo ~ 0.6875 [261, l27|, respectively. This feature of 7oo provides a distinctive signature for 
/(T) as opposed to the ACDM and DGP models. So, one can distinguish the /(T)gravity 
from other cosmological models by the growth index of matter density perturbations. 




FIG. 1: while ^ImO = 0.272, g and 7 are displayed as a function of z for f{T) gravity respectively. 



IV. GROWTH INDEX PARAMETRIZATION AND OBSERVATIONAL CON- 
STRAINTS 



From the right panel of Fig. ([T]), it is easy to see that, for any n, 7 is not a constant, 
especially in the redshifts region {z < 2) where some observational data points are ob- 
tained. So, it is unreliable to discriminate different models with these observational data 
if the growth index is treated as a constant. The growth index 7 should be a function of 
z and we may parameterize it. Since Eq. fl2U]) can be reexpressed as 



3a 



[n 



1)E 



2n-2/ 



mO) 



dg 



1 -nE2«-2(l 
3 O 



mO ) 



+ 9" + 9 



3 1-E 



2n-2/ 



1-a 



mO) 



2 1 - nE^ 



n-2( 



^mo) 



2 1 _ «(i-n^o)E^"~2 ' 

2n-l 

substituting Eq. ^ into the above expression, we obtain an equation of 7(2;) 



;i + z)lnfi„7' + fi;^ + (2-37)--(l 



27) 



E 



2n-2/ 



mO j 



l-nE^^-^l-n^^) 



2l- 



n(l~n^o)g2" 



2n-l 



(22) 



(23) 



Here, we consider a parametrization form of 7(z) proposed in Ref. 41| 



7(2;) = 7o + 71 



(24) 



which gives a very good approximation of 7(z) for the wCDM and DGP models. The 
error is below 0.03% for the ACDM model and 0.18% for the DGP model for all redshifts 
when f2mo = 0.27. In f{T) gravity, for any given values of n and fimo^ we can obtain 
the value of 70 through the value of determined from solving Eq. fl22l) numerically. 
Substituting Eq. (l24l) into Eq. (123!) . we can get the expression of 71 which is a function 
of redshifts z. For simplicity, we take the value of 71 at 2; = 0, 



71 



-1 \-i 



070 



2 1 _ ra(l-Slmo) 
2n-l 



Q X , 3 fi„o(l - 270) 
(2 - 370) + -- 



21 -nil - a 



mO ) 



■ (25) 



Thus, we obtain the possible regions of 70 and 71 for 0.20 < Qmo ^ 0.35. The results are 
shown in Fig. ([2]). We find, for a given value of Qmo, that the values of 70 and 71 in f(T) 
gravity {n 7^ 0) are larger than those in the ACDM model (n = 0), and a larger value 
of n gives larger values of 70 and 71. The reason for these is that the effective gravity 
in f{T) theory is weaker than that in general relative, and a larger n leads to a weaker 
effective gravity. These features provide distinctive signatures for f(T) gravity as opposed 
to the ACDM model. Therefore, in principle, we can discriminate the /(T) model from 
the ACDM model merely through the values of 70 and 71 if one can obtain their accurate 
values from the observation data. 



0.66 r 

0.64 ' 
0.62 



0.20 0.22 0.24 0.26 0.28 0.30 0.32 0.34 



0.20 0.22 0.24 0.26 0.28 0.30 0.32 0.34 



FIG. 2: 7o and 71 are displayed as a function of QmO for n = 0,0.1,0.25 respectively. 
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Before studying the observational constraints on 70 and 71, we need to examine how 
well the fim^'' with 7(2;) taking the parametrization in Eq. flMj) approximates the growth 
factor g. Numerical results are shown in Fig. (|3]) with Qrno = 0.272. From this figure, we 
see that the error is below 0.45% for n = 0.25 and below 0.1% for n = 0.1. So, ilm 
approximates the growth factor g very well both at low and high redshifts in /(T) theory, 
and we can use all data points to constrain this parametrization. 



0.010 




FIG. 3: The relative difference between the growth factor g and Qm with 7 = 70 + lij^ 
and $7mo = 0.272. The dashed, dotted and sohd curves show the results of n = 0.25,0.1,0, 
respectively. 



In order to obtain the observational constraints on 70 and 71, we first need to determine 
the value of flmo and n from the observations. Here we use the results obtained from the 
combination of the latest Union2 Type la Supernova (Sne la) set, the BAO from the 
SDSS data and the Cosmic Microwave Background (CMB) radiation [12]. At the 95% 
confidence level, f2mo = 0.2721q;q32, n = 0.041;q;33 for the power law model. With these 
best fit values and Eqs. (120} |25|) . we obtain that the corresponding theoretical values of 
7o and 71 are 70 = 0.564 and 71 = —0.0123. 

Now we discuss the observational constraints on 70 and 71 from the growth factor 
data. Here 12 data points given in Table [I] are used. Let us note that although the data 

9 



given in Refs. 42|, |43| are measured without 'any' bias, other data points are obtained by 
assuming a flat ACDM model with Q^o taking a specific value, for example, flmo = 0.25 
or 0.30. So, caution must be exercised when using these data. With this caveat in mind, 
it may still be worthwhile to apply the data to fit the models 271, |44, |45| . With the best 
fit values of Qmo and n, we can obtain the constraints on 70 and 71 from the observations 
by using the following equation 



2 _ \ ^ [yobsyZi) — ^Lm \ 
i=l 9« 



(26) 



where agi is the la uncertainty of the g{zi) data. We find that the best fit values are 
7o = 0.809 and 71 = —0.942. The allowed regions at 1 and 2(j confidence levels are shown 
in Fig.(|l]), from which, one can see that the power law model in /(T) gravity is consistent 
with the observations, since the theoretical values of 70 and 71 obtained by using the best 
fit values of VLmo and n for the power law model are in the allowed region at la confidence 
level. 



z 


Qohs 


References 


0.15 


0.49 ±0.1 


[46,42 


0.35 


0.7 ±0.18 


[48] 


0.55 


0.75 ±0.18 


[49] 


0.77 


0.91 ±0.36 


[46,42 


1.4 


0.9 ±0.24 


[50] 


3.0 


1.46 ±0.29 


[51] 


2.125 - 2.72 


0.74 ±0.24 


[42] 


2.2-3 


0.99 ± 1.16 


[43] 


2.4-3.2 


1.13 ± 1.07 


[43] 


2.6 - 3.4 


1.66 ± 1.35 


[43] 


2.8 - 3.6 


1.43 ± 1.34 


[43] 


3-3.8 


1.3 ± 1.5 


[43] 



TABLE I: The summary of the observational data on the growth factor g. 
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ro 

FIG. 4: The la and 2a contours of 70 and 71 by fitting the power law model in f{T) gravity 
to tlic growth rate data. The point denotes the theoretical values of 70 and 71 with the ^mO, n 
taking the best fit values. 



V. CONCLUSION 



In this paper, we study in detail the growth index of matter density perturbations for 
the power law model in f{T) gravity. Using a parametrization 7(2;) = 70 + 7i for the 
growth index 7(2;), which gives a very good approximation of 7(-2), wc find that the value 
of 7o and 71 in /(T) gravity arc larger than those in the ACDM model, and a larger value 
of n gives larger values of 70 and 71. This feature may provide a signature for f{T) gravity 
distinctive from the other models, such as the ACDM and DGP. Finally, we discuss the 
constraints on 70 and 71 from the observational growth factor data and find that, at lo" 
confidence level, the power law model in f{T) gravity is consistent with the observations 
since the theoretical values of 70 and 71 obtained by using the best fit values of and 
n are in the allowed region. 
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